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ABSTRACT 

Precisely measured astrometric parameters are integral to successful pulsar timing campaigns. They 
are commonly measured by fitting the astrometric parameters of a deterministic timing model to a 
series of pulse times of arrival (TOAs). TOAs measured to microsecond precision over several- year 
spans can in this way provide astrometric parameters precise to sub-milliarcsecond levels. However, 
pulsars do not pulsate in a deterministic fashion. Many display significant amounts of red spin noise. 
Furthermore, a stochastic background of gravitational waves can lead to red noise-like structure in 
TOAs. We investigate how noise of different spectral types is absorbed by timing models and leads to 
significant estimation errors in the astrometric parameters. Independent of timing, very long baseline 
interferometry (VLBI) is capable of providing sub-milliarcsecond astrometric parameters for pulsars. 
We find that incorporating VLBI astrometric measurements into the timing models of pulsars for which 
only a couple of years of timing data exist will lead to more realistic assessments of red spin noise, 
yield more accurate astrometric parameters, and could enhance the amplitude of certain gravitational 
wave signatures in post-fit timing residuals by factors of 20 or more. 
Subject headings: astrometry - gravitational waves - pulsars: general - reference systems 



1. INTRODUCTION 

Successfully timing a pulsar relies on the extraordinar- 
ily stable evolution of the rotational phase of that pulsar 
and the fact that its phase can be tracked as it beams 
pulses of light towards Earth at nearly fixed points in 
its phase. The pulses traveling towards Earth leave the 
pulsar at intervals indicating the phase evolution of the 
pulsar, but those intervals are modulated before reach- 
ing Earth by many astrophysical processes. To regain 
the phase information of the pulsar, one must develop a 
detailed mathematical model (a timing model) that ac- 
counts for all the deterministic physical effects that cause 
pulses to arrive at Earth at intervals differing from those 
with which they were emitted. The main data prod- 
ucts of a pulsar timing campaign are a set of parameters 
describing the astrophysical phenomena being modeled 
and a series of timing residuals — the differences between 
measured pulse arrival times and the predictions of the 
timing model. A successful timing model paired with 
a pulsar rotating in a predictable w ay will yie ld timing 
residuals consistent with white noise ([Manchest er & Tay- 
lor|1977[pandford et al.|1984[[Lorimer fc Kramer|2005( )r 



'he precise astrometric parameters ol a pulsar (right 
ascension and declination (a, i5) = 0, proper motion 
{^a, ^ls) = and parallax tt) are essential to success- 
ful timing models. An error in position of just a tenth of 
an arcsecond can lead to annual oscillations in residuals 
of over 100 /iS in an iplitude (see Figure 1). With soft- 



ware like TEMP02 (Hobbs et al. 2006), timing models 



are refined through least squares iitting to minimize the 
amplitude of such signatures in the residuals. For a par- 
ticularly stable subset of recycled (millisecond) pulsars 
(MSPs), measuring pulse times of arrival (TOAs) with 
microsecond precision or better and refining the timing 
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model parameters to minimize the RMS o f the residuals 
can yield sub-m illiarcse cond astrometry ( Splaver et al. 
2005||Hotan et al..2006,|Verbiest et al.|2008p . Precise as- 
trometric measurements are also useful m studies of the 



Galactic neutron star (NS) population ( Cordes fc Cher 



noffil9971 [Toscano et aL][TM9l piucher-Giguere Kasp: 
2007[), NS e lectromagnetic emission mechanisms fDeller 
et al.|200 9b |, the Gala ctic distribution of interstellar gas 
( Corcles & Lazio|2002 1 , and the evolution of compact bi- 



naries containing pulsars, specifically to correct for the 



Shklovskii effect (Shklovskii 



19691 IDeller erar][2008l ). 



As pulsars are point-like radio sources, their astromet- 
ric parameters can also be obtained with very long base- 
line interferometry (VLBI). VLBI is currently capable 
of sub-milliarcsecond astrometric measurements of pul- 
sars on par with the best timing-based results. Addi- 
tionally, VLBI astrometry can arrive at precisions re- 
quiring approximate ly five years of timing data in less 
than two years (e.g., Chatterjee et al.|2009 ). To enhance 
the observation signal-to- noise ratio (s/N), pulsar gating 
(gathering signal only when the pulsar is beamed towards 
Earth) is used. The Very Long Baseline Array (VLBA) 
allows for full-time VLBI and better-controlled system- 
atics compared to other VLBI efforts owing to the identi- 
cal antennas in the array. With gating, the VLBA, and a 
growing population of sources useful for calibration, the 
astrometric precisions achievable with VLBI continue to 
improve. 

Canonical pulsars (CPs) and, to a lesser extent, MSPs 
display what is known as red spin noise. Red spin noise is 
likely caused by complex dynamics at the crust-core in- 
terfaces or in the magnetospheres of pulsars. It manifests 
itself as noise in timing residuals with an approximately 
power-law power spectrum P{f) cx (more power at 
lower frequencies, hence the moniker "red") and can be 
characterized by stochastic shifts in the period or period 
derivative of a pulsar. The stochasticity of red spin noise 
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Fig. 1. — The characteristic signature in timing residuals caused 
by errors in position, proper motion, and parallax (from top 
to bottom respectively). For a pulsar with the coordinates of 
PSR 1713-1-0747, these are the residuals associated with 0.1 arcsec- 
ond errors in 5 and vr and 0.1 arcsecond per year errors in fig. The 
oscillations in the top two plots have a 1 year period and the oscil- 
lations in the bottom plot have a 6 month period. The amplitudes 
of these curves scale linearly with the errors in the astrometric 
parameters. 



makes it impossible to incorporate into a deterministic 
timing model. In contrast to white noise, the excess of 
low frequency power in red spin noise introduces corre- 
lations over finite timescales. These correlations give red 
spin noise a non-trivial shape that can vary substantially 
from realization to realization. 

Along with red spin noise, gravitational waves (GWs) 
induce structure in timing residuals. There are cur- 
rent efforts t o detect the signature of isolate d, contin- 
uous so urces ( Lee e t al. 2011; Ellis et al. 2012) and burst 
sources (Finn"SF Tommcn 2010 
Haasteren & Levin 2010 



Pshirkov et al.|2010 
Cordcs & Jenct 201 



van 
of GWs 



in timing residuals, but among the most promising types 
of GW signals from an initial detection perspective is 
that of a stochastic background (SB) created by an en- 
semble of merging supermassive black holes (SMBHs) in 
the distant universe. The SB is characterized by strains 
with an amplitude spectrum hc{f) (x f~^^^ that mani- 



spectrum P(f) oc / ^^/^ ( 


Jenet et al.|2006 


Sesana et al. 


2008 


Cordes & Shannon|2U12 ). Since red spin noise and 



indices and could be of comparable magnitudes, it is ex- 
ceedingly difficult to make a detection of this typ e of GW 
in the res i duals of a single pulsar. Fortunately, |Hellings| 



& Downs ( 1983 1 predict that a specific pattern of corre 
lations unique to GWs can be observed in the residuals of 
an array of pulsars in the presence of a S B of GWs . The 



European P ulsar Timing Array (EPTA; van Haasteren 
et al ||2011|), the P arkes Pulsar Timing Array (PP'i'A; 



Yardley et al.|2011 ), and the North American Nanohertz 
Observatory for U ravitational Waves (NANOGrav; 'De-J 



morest et al.|2012 ) are timing arrays of pulsars and look- 
mg tor, among other things, the type of correlations in 
their residuals predicted by Hellings and Downs. 

If one is trying to study the red spin noise of a pulsar 
or is trying to detect GWs with a PTA, the model fitting 
applied in timing studies may undermine those efforts. 



From a statistical standpoint, every fit parameter reduces 
the degrees of freedom in the data set. Furthermore, if 
one produces an initial timing model with approximately 
correct parameters, the timing residuals will contain as- 
trophysical information, but will also contain the signa- 
tures of errors in the timing model. Errors in each of the 
timing model parameters leave characteristic signatures 
in the timing residuals that have associated structures in 
frequency space. Minimizing the amplitude of the signa- 
ture of a timing model parameter error indiscriminately 
removes power from the residuals at the freq uencies as 
sociated w ith that particular signature (e.g., Blandford 
eFaLlflOSil ). The power that is removed from the pre-fit 



residuals is absorbed by the timing model and causes the 
parameters of the timi ng model to shif t, possibly away 



from their true values ( Ellis et alT][2011 1 . 

Because the spin noise intrinsic to many pulsars and a 
SB of GWs introduce red processes into timing residuals, 
we investigate the errors induced in the astrometric pa- 
rameters of a timing model when they are fit to a series 
of residuals containing red noise power. In so doing, we 
demonstrate how red processes can be partially absorbed 
into a timing model and thus attenuated in the residuals. 
The deleterious effects of model fitting that we demon- 
strate can be side-stepped if the timing-independent as- 
trometry of VLBI is incorporated into future pulsar tim- 
ing models. 

In Section 2, we discuss astrometry as it pertains to 
pulsar timing models, and in particular, derive the sig- 
natures in timing residuals associated with errors in po- 
sition, proper motion, and parallax. In Section 3, we 
describe simulations we conducted to assess the effects 
of fitting for astrometry in the presence of red processes 
on the astrometric parameters of a timing model and on 
the signature of the red process itself. We follow this 
with a discussion of our results in Section 4. In Section 
5, we discuss the real- world feasibility and benefits of 
incorporating VLBI astrometry into pulsar timing cam- 
paigns. In Section 6, we discuss several applications for 
VLBI astrometry in pulsar science that are not directly 
tied to gravitational wave detection. Finally, in Section 
7, we summarize our results and provide some concluding 
remarks. 

2. ASTROMETRIC TERMS IN TIMING MODELS 

In this section, we discuss the roles of position, proper 
motion, and parallax in the development of timing mod- 
els. For a more thorough descripti on of this m aterial, see 
[Backer fc Hellings| (11986) and |Edwards et al.| (|2006). The 
pulsar position anoproper motion principally affect the 
Roemer delay. The Roemer delay is the delay between 
the arrival of a pulse at an observatory and the arrival of 
that pulse at the solar system barycenter (SSB). It can 
be written as An = —r • Rp/c where r is the vector 
pointing from the SSB to the geocenter (we consider er- 
rors in the topocentric to geocentric transformation to be 
negligible), and Rp is the unit vector pointing from the 
observatory to the pulsar at the time of the observation. 

We use equatorial coordinates a and 6 (referenced to 
the geocenter) in all other sections of this paper, but here, 
we use an ecliptic coordinate system (centered at the 
SSB) with longitude A and a latitude /3. Let A increase 
with the orbital motion of Earth from an x-axis that 
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points from the Earth to the SSB at the vernal equinox. 
Place the z-axis such that the Earth orbits in a right- 
handed fashion. Exploiting the relation r + Rp = rp, 
where rp is the vector pointing from the SSB to the pul- 
sar, and approximating the orbit of the Earth as circular, 
we write 

r = —xRcos (ut) — yRfim{ujt) (1) 
rp = xD cos /? cos A -t- yD cos (3 sin A + zD sin /3 (2) 
Rp — x[D cos P cos A -f i? cos (wt)] + 

y [D cos P sin A -t- i? sin (cjt)] + zD sin /3, (3) 

where i? is 1 A.U. and uj is 27r/(l yr). If the position and 
proper motion of the pulsar are determined at i = and 
if any acceleration of the pulsar (from motion in a binary. 
Galactic rotation, etc.) is ignored, the position can be 
written as j5{t) = /3o + A'/gi and \{t) = Aq + jixt where 
jjLx and /i^ are the components of the proper motion. 
By expanding to linear order in the small parameters 
fix, ^p, and n = R/D and dropping any constant terms, 
the Roemer delay can be written as 



[cos/3o cos {ujt — Aq) 

— fipt sin Pq cos {tut — Ao) 

— fixt cos /?o sin (ujt — Aq) 

— — cos^ Po cos {2ujt — 2Ao 



(4) 



Errors in the astrometric parameters, assumed small, 
produce timing variations that are sinusoids with peri- 
ods of one year and six months (position and parallax, 
respectively) and a linearly growing sinusoid with a one- 
year period (proper motion). This is illustrated in Figure 
1. The amplitudes of these oscillations scale linearly with 
the magnitude of the astrometric error. If a set of residu- 
als display structure resembling a superposition of these 
signatures, fitting the model to the residuals will lead to 
adjustments of the astrometric parameters of the model. 

3. SIMULATIONS OF NOISE-INDUCED ASTROMETRIC 
ERRORS & RED NOISE ABSORPTION 

We simulated three sets of noiseless arrival times for 
a pulsar with the astrometry, spin frequency, and spin- 
down rate of PSR J1713-1-0 747. We did so w ith the 
"fake" plug-in for TEMP02 ( |Hobbs et al] [2006)). In the 
first set of pulse arrival times (Set A) we simulated one 
TOA every 3 days over a 500 day span while in the sec- 
ond set (Set B) we simulated one TOA every 21 days over 
a 500 day span. In the final set (Set C) we simulated one 
TOA every 21 days over a 3,000 day span. 

We added noise of different amplitudes and spectral 
characteristics to the simulated TOAs. We considered 
red noise with a power spectrum P{f) oc f^'^ with spec- 
tral indices 7 = 0, 3, and 5. For each 7, we chose 16 
amplitudes for the pre-fit RMS of the noise logarithmi- 
cally spaced between 10 ns and 10 ^s. Once combined 
with a particular realization of noise of a specific spectral 
type and RMS amplitude, we reprocessed the simulated 
pulse arrival times with TEMP02, fitting for the phase 
(0), spin period (P), spin-down rate (P), and astrome- 
try. We recorded the RMS of the post-fit residuals along 
with the astrometric parameters. For each spectral type 



of noise and each pre-fit noise amplitude, we did this with 
100 different realizations of noise. 

We calculated the post-fit astrometric errors as fol- 
lows: 



I = [(a - apost)'^ cos^ S+{S- Spostf 



,1/2 



(5) 



5fl = [{fla - fia,postf COS^ 5 + {fis - fiS,postf] (6) 

(57r = Itt - TTpostl. (7) 

Here, a, (5, /Iq, fis, and tt are the astrometric parameters 
used to simulate the initial TOAs; the quantities with a 
subscript "posi" are the post-fit estimates for them. 

4. DISCUSSION OF SIMULATION RESULTS 

In Figure 2, we show histograms of the factors by which 
the pre-fit RMS amplitude of the noise are reduced by 
fitting. For simulation Sets A and B, which correspond 
to 500 days of observations, significant amounts of red 
noise are absorbed into the model. In both cases, the 
redder the noise, the more severe the absorption. In the 
reddest case, the pre-fit RMS amplitude of the red noise 
is reduced by a factor of about 20 in the median case. 
The problem is less severe in the simulations with longer 
observation spans (Set C). Even so, in the reddest case, 
the pre-fit levels of red noise can be reduced by up to a 
factor of five. 

Figures 3, 4, and 5 (corresponding to simulations Sets 
A, B, and C, respectively) show errors in the astrometric 
parameters of the timing model when they are fit in the 
presence of noise. In all cases, the noisier the pre-fit resid- 
uals, the noisier the post- fit residuals and the bigger the 
errors in the astrometric parameters. These figures also 
depict the reduction of the RMS noise level that is caused 
by fitting. Like in Figure 2, it is apparent that the ab- 
sorption of red noise is a more severe problem for the 500 
day simulations in Sets A and B than it is for Set C; this 
can be straightforwardly explained. Fitting for phase, 
period, and period-derivative removes a quadratic from 
the pre-fit residuals. The post-fit residuals will thus con- 
sist of cubic and higher-order contributions. The higher 
order contributions will be suppressed in red processes 
that have less power in higher frequencies. A zero-mean 
polynomial dominated by a cubic contribution can read- 
ily mimic a sinusoid with a one year period if the data 
set has a length of approximately one year. If the data 
set is significantly longer than one year, this is not such 
a problem. Similarly, if one has only 6 months of data, 
quadratic subtraction may leave a shape in the residuals 
resembling the 6-month-period sinusoid associated with 
errors in parallax. 

In Figures 3, 4, and 5, there is a stark difference be- 
tween 7 = (white noise) and the 7 = 3 and 5 cases. 
With white noise, the 100 points associated with each of 
our 16 levels of input pre-fit RMS noise fall into signif- 
icantly more compact regions of the plane for position, 
proper motion, and parallax. These regions are partic- 
ularly compact in the horizontal dimension. All realiza- 
tions of white noise are comparably non-covariant with 
the signatures of astrometric errors in the timing model. 
On the other hand, some realizations of red noise can 
be more covariant with the signatures of astrometric er- 
rors in the timing model than other realizations. Some 
realizations of red noise can be more readily attenuated 
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Fig. 2. — Ratios of post-fit RMS to pre-fit RMS for simulated ar- 
rival times. Each of the three histograms in each panel corresponds 
to input noise of a different spectral index (7 = 0,3,5). Each his- 
togram contains 1600 counts for the 100 realizations of noise at 
each of 16 pre-fit rms amplitudes. Redder noise is more readily 
absorbed into the astrometric parameters of a timing model, but 
the amount of absorption is reduced as observation spans grow 
significantly longer than 1 yr. 




10' 10' 10" 

Noise RMS \fis] (7 = .5} 



Fig. 3. — Astrometric errors (offset between true astrometric pa- 
rameters and post-fit estimates of astrometric parameters) versus 
post-fit RMS for simulation Set A (500 days of TOAs taken at 3 
day intervals). The three panels of the three plots each contains 
1600 points corresponding to the 100 realizations of noise at 16 
input levels of pre-fit noise. Moving downwards through the three 
plots, the spectral index of the additive noise increases from 0, to 3, 
to 5. The yellow points correspond to an input pre-fit rms of 0.01 
fis (this is indicated by the triangular magenta indicator above 0.01 
fis). Similarly, cyan, magenta, and blue points correspond to input 
pre-fit rms values of 0.1, 1, and 10 /is repsectively. The top, middle, 
and bottom panel of each plot indicates the noise-induced post-fit 
error in the position, proper motion, and parallax respectively. 




Fig. 4. — As in Figure 3, but for simulation Set B (500 days of FiG. 5. — As in Figure 3, but for simulation Sot C (3000 days 

TOAs taken at 21 day intervals). TOAs taken at 21 day intervals). 
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in post-fit residuals than others by modifying the astro- 
metric parameters of the timing model; this leads to the 
increased horizontal scatter in simulations involving red 
noise. 

5. INCORPORATING VLBI ASTROMETRY INTO PTAS 

Given a sufficient S/N and sufficient calibrators, mod- 
ern VLBI can provide absolute position estimates relative 
to the International Celestial Reference Frame (ICRF) 
to 1 mas precision after 1 year, proper motion estimates 
to 0.01 mas yr~^ precision after 2 years, and parallax to 



from pulsar to pulsar, the errors in the frame transfor- 



mation can be reduced by a factor cx N 



-1/2 



At the 



0.01 m as precision after three years (e.g., Chatterjee et al. 
20091. From Figures 3, 4, and 5, we can infer the post- 



tit KIVIS noise amplitudes at which noise-induced errors 
in the astrometric parameters reach levels comparable to 
the precisions achievable with VLBI. A summary of our 
conclusions is presented in Table 1. The best-timed pul- 
sars to date have post-fit RMS residuals of several tens 
of nanoseconds over several yea r spans ( ^ Yardlcy ct_al^ 
'2011 van Haasteren et al. 2011 [Demorest et af. |26l2 ). 
In all cases in I'able 1 (except tor position estimates m 
our 3000 day simulated data sets), VLBI can provide 
better astrometric measurments than timing techniques 
for nearly all known pulsars, even the best-timed ones. 

5.1. Reference Frame Considerations 

Astrometric me asurements with VLBI are made rela- 
tive to the ICRF (Ma et al."1998'), a geocentric, quasi- 
inertial frame referenced to distant quasars. Pulsar tim- 
ing is carried out in a dynamic reference frame cen- 
tered on the SSB and oriented relative to the orbital 
plane of Earth; this coordinate system is set by a so- 
lar system cphcmcris such as the Jet Propulsion Labora- 
tory's DE405 (Standish 1998). Position measurements of 
some pulsars with VLBI and pulsar timing are inconsis- 
tent, indicating that the quoted measurement errors are 
too small. To demonstrate the inconsistency between 
VLBI and timing-based astrometry, in Table 2, we show 
VLBI and timing-based position measurements of pulsars 
J0437-4715 and J1713-t-0747. The two timing-based po- 
sition estimates for pulsar J1713-I-0747 in Table 2 are 
also inconsistent with one another. This inconsistency 
is also d ue, at least in part, to reference fra me differ 
ences as 



Splaver et al. (20051 use t he DE200 dStandishl 

19901 solar system ephemeris while Hotan et aT ( 2006[ ) 



use the more accurate DE405 ephemeris. Proper motion 
and parallax measurements are less susceptible to such 
frame-related offsets as they depend more on the changes 
in celestial coordinates over time than the absolute co- 
ordinates. However, as Eq. 4 demonstrates, timing cor- 
rections associated with proper motion and parallax de- 
pend on the position of the pulsar, so reference frame 
mismatch will influence the degree to which VLBI-based 
and timing-based proper motion and parallax measure- 
ments can be made to agree. 

By measuring the absolute position of a pulsar with 
both timing and VLBI at approximately contemporane- 
ous epochs, it is possible to constrain the two-angle ro- 
tation needed to transform between the ICRF and the 
relevant pulsar timing r eference frame to m illiarcsecond 
precision or better (e.g., [Bartel et al.||1996[ ). If compara- 
bly precise VLBI and timing astrometry is measured for 
Np pulsars and the measurement errors are uncorrelated 



outset of a timing campaign on a sample of newly dis- 
covered pulsars, pre-existing measurements such as those 
presented in Table 2 can be used to refine the frame tie 
to an extent that would allow VLBI astrometry to im- 
mediately be used in timing models in place of estimates 
derived from timing data. With lar ge VLBI a strometry 
campaigns like the PSRtt project (De ller et al.|2011 ) cur- 
rently under way, it is likely that all pulsars m P'l As will 
soon have VLBI astrometry available. As TOAs are col- 
lected, timing models should be refined using standard 
techniques, but the astrometric parameters of the model 
will not have to be fitted. 

Once the time baseline has significantly exceeded one 
year, the PTA data can be made independent of the pre- 
existing frame tie by combining the timing data and the 
VLBI astrometric measurements and doing one simulta- 
neous, global fit over all pulsars and the frame tie. Us- 
ing timing residuals from all pulsars in a PTA to con- 
strain the frame-tie parameters will require incorporat- 
ing currently existing tools like TEMP02 into a more 
elaborate framework that can simultaneously fit the tim- 
ing models of all the pulsars in the array as well as the 
frame-tie parameters with appropriate weighting to ac- 
count for different uncertainties. For a PTA with 17 pul- 

fitted 1382 parameters. 



sars, 



Demorest et al 



(2012) fitted 1382 parameters. Of 
these parameters, llYldealt with a time- varying disper- 
sion measure or pulse profile evolution across the band. 
With more thorough multi-frequency observation, many 
of these parameters will become unnecessary. The re- 
maining frequency-dependent parameters will be deter- 
mined by features unique to the residuals of individual 
pulsars and will prove highly non-covariant with param- 
eters that affect the quality of the global fit. Of the re- 
maining 211 parameters, 85 are astrometric parameters 
that could be provided by VLBI if 2 additional frame-tie 
parameters are included in the fitting procedure. 



5.2. Incorporating New Pulsars Into PTAs While 
Utilizing VLBI Astrometry 

Suppose a new MSP is discovered and its timing stabil- 
ity and brightness make it a prime candidate for PTA sci- 
ence. In the first couple of years of observations, if timing 
and VLBI studies begin concurrently and are combined 
appropriately, attenuation of the amplitude of red pro- 
cesses in the timing residuals caused by model fitting can 
be dramatically reduced. T o quantify this, w e introduce 
some formalism taken from B landlord et al.| (jT984). The 
timing residuals, R{t), are fit to a linear combination of 
basis functions associated with the parameters being fit 
for. To simplify matters, consider only eight parameters 
as follows: 



R{t)^ 


5I^'^V'^(0 

1=1 




Mt) = 1 




cos {ut) 






t cos {ut) 




^7{t) = 


t sin (iot) 


ipiit) = sin (ujt) 


Mt) - 


sin {2ujt). 



(8) 



7 



TABLE 1 



NDays 


gap 


NTOA 


7 


Ag{l mas) 


Ap(0.01 mas yr"!) 


A^(0.01 mas) 


(days) 


(days) 






(ns) 


[ns) 


(ns) 


500 


3 


169 


5 


< 1 


<1 


<1 


500 


3 


169 


3 


< 2 


<1 


<1 


500 


3 


169 





35 


<10 


<10 


500 


21 


25 


5 


1 


< 1 


< 1 


500 


21 


25 


3 


2 


< 1 


<1 


500 


21 


25 





8 


<10 


<10 


3000 


21 


145 


5 


100 


5 


2 


3000 


21 


145 


3 


120 


7 


3 


3000 


21 


145 





200 


< 10 


<10 


Note. 


— Levels of post-fit 


noise RMS indicative of astrometric errors on par with 



precisions achievable with VLBI astrometry. Moving from left to right, the first two 
columns indicate the number of days worth of TOAs and the gap (in days) between 
TOAs for the set of simulated TOAs used for the simulations associated with that 
row. The third column is the total number of simulated TOAs in that set. The fourth 
row indicates the spectral index of the power spectrum of the noise added to the 
pristine set of TOAs. The last three columns indicate the RMS amplitude of noise 
in the post-fit residuals at which the noise-induced error in position, proper motion, 
or parallax is big enough that VLBI would provide more accurate astrometry. If 
a "less-than" sign is used, the interpretation is that post-fit RMS noise well below 
that level is indicative of errors bigger than precisions available from VLBI. 



TABLE 2 



PSR 


Measurement Type 


Aa 


AS 


Reference 




(s) 


(") 





J0437-4715 



VLBI 
Timing 



0.883250(3) 
0.883185(6) 



0.031863(37) 
0.034033(70) 



Deller et al.| 2008| 



Verbiest et al.M^U08 



lilqolL 



J1713+0747 


VLBI 


0.5306(1) 


0.519(2) 


Chatterjee et al. ( 2009|l 




Timing 1 


0.5307826(7) 


0.52339(3) 


spiaver et al.^ ^::!U05|I 




Timing 2 


0.53077(1) 


0.5228(2) 


Hotan et al.|(2uub|l 



Note. — Astrometric positions measured with VLBI and pulsar timing for PSR 
J0437-4715 and PSR J1713-f0747. For PSR J0437-4715, Aa is the offset from right as- 
cension 04''37'"15^ and AS is the offset from declination -47°15'09" at MJD epoch 54100. 
For PSR J1713-I-0747, Aa is the offset from right ascension 17'' 13^49= and AS is the offset 
from declination 07°47'37" at MJD epoch 52275. 



The basis functions i^iit) are associated with 
(f), P, P,a, fia, fj^s, and tt respectively. The coeffi- 
cients Xi are functions of the true spin and astrometric 
parameters of the pulsar and the fractional errors in the 
timing model estimates of them. Assuming N years of 
observations have been made, the basis functions can 
be converted to an orthonormal set of basis functions, 
ip'iit), using the Gram-Schmidt process with the inner 
product 



N/2 
'N/2 



(9) 



We then compute the Fourier transform of each basis 
function, and combine them into a transmission 

function, T{f) as follows: 

fN/2 

i'Kf)^ V:Wexp(27ri/t)di (10) 

J-N/2 

nf)=i-l^jz^[{m*{f)- (11) 
1=1 

If the structure in the pre-fit timing residuals is due en- 
tirely to some process with a stationary power spectrum 



P{f), the mean-square post-fit residual coming from fre- 
quencies between and fn is given by 



i?2: 



fn 



fL 



T{f)PU)df. 



(12) 



The transmission function is a filter on the pre-fit vari- 
ance and dictates what fraction of power in the pre-fit 
residuals at a particular frequency makes it through the 
fitting procedure and into the post-fit residuals. 

If sufficiently precise VLBI astrometry is available and 
reference frame mismatches have been dealt with, cer- 
tain parameters do not need to be fit for. The process 
described above can be repeated with any subset of the 
basis functions listed in Eq. 8. We assume that after the 
first year of timing and VLBI observations, VLBI esti- 
mates for position are precise enough that they do not 
need to be fit for. After 2 years, the VLBI estimates 
for proper motion reach sufficient precision, and after 3 
years, VLBI can determine parallax well enough that no 
astrometric fitting is necessary. In Figure 6, we demon- 
strate how the transmission function varies in the first 5 
years of observation if VLBI astrometry is or is not in- 
corporated into timing efforts according to this schedule. 

The sensitivity enhancements when VLBI astrometry 
is used are apparent in Figure 6. Using VLBI astrome- 
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Fig. 6. — Transmission functions for 1, 2, 3, and 5 year data sets 
with and without using astrometric measurements from VLBI. The 
legends indicate the spin and astrometric parameters being fit to 
produce the associated transmission function. In the 5 year trans- 
mission function in which all 8 parameters are fit, the prominent 
absorption feature at the origin (/ = yr~^) is associated with the 
quadratic fit to the spin parameters of the pulsar, the feature at 1 
yr~^ is associated with fitting out the oscillations associated with 
errors in position and proper motion, and the feature at 2 yr~^ is 
associated with fitting out parallax errors. 







— 1 .yr 

— 7yt 

— 3}T 
1} JT 




"'012345 



Fig. 7. — The ratio of post-fit residual variance using VLBI as- 
trometry to that without using VLBI astrometry versus the spec- 
tral index of a red process in the TOAs. Each curve is associated 
with one panel from Figure 6, as is indicated by the legend. It is 
assumed that astrometric parameters from VLBI are incorporated 
into the timing model according to the same schedule used to make 
Figure 6; after 1 year, VLBI provides 0, after 2 years, it also pro- 
vides 11, and at the 3 and 5 year marks, VLBI is providing all 5 
astrometric parameters. The quantity being plotted is Ry from 
Eq. 13. 



try, it will be possible to detect coherent sources of GWs 
with frequencies of approximately 1 yr^^, a prospect that 
is entirely hopeless without VLBI. Aside from this very 
specific type of source, because the bulk of the sensitiv- 
ity enhancement occurs at low frequencies (below about 
3 yr~-^), the sensitivity of PTAs to red processes could be 
greatly enhanced. To quantify the enhancement to red 
process sensitivity, for a range of spectral indices 7, we 
compute the following variance ratio: 



//" TyLBlU)ndf 



Rvh)-- 



ll"T{f)f-^df 



(13) 



If, 



where TvLBiif) is the transmission function when VLBI 
astrometry is used and T(/) is the transmission function 
when all the astrometry is fit for. We set the low fre- 
quency cutoff /l = 0. The low-frequency absorption fea- 
ture common to all the transmission functions in Figure 6 
approaches zero as for small /, so each of the integrals 
in Eq. 13 will converge for 7 < 7. We set the high fre- 
quency cutoff fn = 9 yr~^ corresponding to the Nyquist 
frequency if 18 equispaced measurements are made each 
year. The results of this calculation are shown in Figure 
7. As 7 increases (steeper power spectrum), the increased 
amplitude of the transmission function at low frequencies 
when VLBI information is used becomes increasingly im- 
portant so as to not filter out substantial amounts of the 
power through fitting. The enhanced sensitivity is most 
pronounced after two years with red processes having a 
spectral index very near 13/3, the value expected from 
a SB of GWs. The enhancement can be greater than a 
factor of 20 in this case. This is consistent with our find- 
ings from Figure 2 that for data sets 500 days in length, 
red noise with a spectral index of 5 is attenuated by fit- 



ting by a factor of 20 in the median case. Once the data 
set becomes longer than three years, the transmission 
functions with and without VLBI astrometry differ only 
by the prominent absorption features at 1 and 2 yr~^. 
The width of these features A/ is approximately propor- 
tional to f /N so the feature at / = 1 yr~^ has a width 
A/ « N^^ yr^^ and the feature at / = 2 yr^^ has a 
width A/ w 2A^^i yr-i. 

6. ADDITIONAL VLBI AND TIMING SYNERGIES 

So far, we have emphasized the significance VLBI as- 
trometry can have in PTA efforts to detect GWs. Here, 
we briefly discuss additional types of pulsar timing sci- 
ence that could be aided by VLBI astrometry. 

6.1. Relativistic Evolution of NS-NS Binaries 

Pulsars in compact binaries show post-Keplerian (PK) 
evolution described by parameters that include the rate 
of periastron advance, the gravitational redshift, and the 
rate of orbital period decay. Within any theory of grav- 
ity, these PK parameters are functions of the unknown 
masses of the two objects in the system and the mea- 
surable Keplerian parameters. Measurements of any two 
independent PK parameters yield estimates of the two 
masses. Measurements of addition al PK param eters pro- 
vide tests of that 
of the binary syste: 



theory of gravity (| Stairs | 2003|| . Studies 
terns PSR B1913-(-16 jHulse & Taylor 



19751 IWeisberg et a l.|[20Tol) and PSRs 

— ■ \2m\ E — ^ ' • 



Taylo: 
U737-3U39A/B' 

( [Kramer et ai.|200tjp k ave m this way provided stringent 
tests of general relativity (GR) and, through measure- 
ments of the rate of orbital period decay, Pb, the best 
indirect evidence for gravitational radiation to date. 

Measurements of Pb are complicated by the fact that, 
in the timing model, a decaying orbital period is co- 
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variant with the relative acceleration of the centers of 
mass of the compact binary and the SSB. Additionally, 
the Shklovskii effect is an apparent relative acceleration 
caused by a changing line of sight to the pulsar that 
depends on th e proper motio n and distance to the pul- 



sar as n^D (Shklovskii] |1969[) . For some binaries, like 
~ kiov; ~ 



inant bias p reventing an accurate measurement of Pb 
(Stairs 2003). If the proper motion and parallax of a 
pulsar are provided by VLBI, the bias can be removed 
peller et al]|200 9a|). 

i<'or B1913+16, despite the existence of over 30 years 
of timing data, VLBI estimates of the astrometric pa- 
rameters will be of higher precision than those derived 
from timing. Pulsars in NS-NS binaries have larger mag- 
netic fields than MSPs a nd therefore are expected to have 
greater spin noise (e.g.. Shannon & Cordes 2010). Over 
approximately 30 years, B1913-I-16 has timmg residu- 
als with an RMS of '--^lOO fis and there was a glitch 
in the rotational frequency of the pulsar in 2003 with 
/S.f / f K, 3.7 X 10"^"'^ (the glitch is indicative of red spin 
noise in B1913+16 at some le v el). D espite decades of 



observations, Weisberg et al. (20101 derive a timing- 



based absolute position measurement for B1913-I-16 with 
'^l mas error. Figure 5 shows that even with over eight 
years of timing data, residuals with an RMS of just 10 /xs 
are indicative of substantial errors in timing-based esti- 
mates of the astrometric parameters (10 mas errors in 
absolute position). Position measurements from VLBI 
that can reach 1 mas precisions in just one year and are 
independent of the timing stability of a pulsar provide 
an improvement over the timing-based position estimates 
for pulsars in NS-NS binaries. For example, Dclle r et al.| 
(l2'009a) combine VLBI astromctry from the Australian 
iTong Baseline Array and timing data for the double pul- 
sar (J0737-3039A/B) and find that in the future, with 
a decade of timing data, measurements of Pt for this 
system will allow a test of GR to 0.01%. 

6.2. Multi-frequency Studies of Newly 
Discovered j-ray MSPs 

Recent years have seen an explosion in the population 
of known MSPs. This is due, in no small part, to the 
Large Area Teles cope (LAT) on t he Fermi Gamma-ray 
Space Telescope. Ray et al. (2012) describe radio follow- 
ups of a number of LAI' pulsar candidates. Of 300 candi- 
dates that were 7-ray sources having spectra character- 
istic of pulsars and signs of variability but no definitive 
source classification, 47 new pulsars were discovered, 43 
of which were MSPs. In one case, by conducting a ra- 
dio follow-up with the Nangay Radio Telescope on an 
unidentified LAT 7-ray source, Guillemot et al. (20121 
discovered the MSP J2043-I-1711. Using an ephemeris 
based on timing data from Arecibo, Nangay, and West- 
erbork, the 7-ray data were folded and pulsations were 
detected in the 7-ray band. Additional follow-up obser- 
vations with Suzaku and the Swift X-ray Telescope failed 
to identify any X-ray counterpart to J2043-I-1711. 

Given the rapid rate of discovery of MSPs by instru- 
ments like LAT, timing observations are limited by the 
availability of telescope time. VLBI astrometry will ac- 
celerate and improve multi-wavelength studies of objects 
like J2043-I-1711. The 7-ray sources from LAT are lo- 



calized to about 1°. If a telescope like the Green Bank 
Telescope can detect a radio pulsar counterpart to a 7- 
ray source, the localization will improve to a few arcmin- 
utes. At this point, standard pulsar timing techniques 
can be used to develop a timing model that will allow for, 
e.g., the detection of pu lsations in 7-rays for J2043-I-1711 
( Guillemot et al.||2012 |. But, until several years of tim- 
ing data are collected, timing will be unable to constrain 
the position of the pulsar sufficiently to allow successful 
identification of a unique X-ray counterpart. Further- 
more, there is a chance that a newly discovered pulsar 
may exhibit high levels of timing noise that substan- 
tially skews timing-based astrometric estimates. After 
the pulsar is localized to within a few arcminutes by an 
initial radio detection, we can instead use interferometry 
(e.g. with the Very Large Array and then VLBI) to con- 
strain the position to milliarcsecond precisions. Efforts 
to detect the pulsar in X-rays can be better guided and 
unique X-ray counterparts can be picked out. Rapid, 
high-precision ast rometry from VLBI wil l facilitate ef- 



forts like that of Guillemot et al. (20121 and will aid 
studies of the pulsar across the entire electromagnetic 
spectrum. 

7. CONCLUSIONS 

Standard model fitting procedures used in pulsar tim- 
ing can falter in the presence of red processes. Whether 
from red spin noise or the activity of a SB of GWs, fitting 
the timing model by minimizing the RMS of the resid- 
uals causes the power in the red process to be absorbed 
into the timing model at the cost of shifting the timing 
model parameters to erroneous values. 

VLBI provides an alternative means of measuring the 
astrometric parameters of a pulsar with sufficiently high 
precision that these measurements can be adopted by 
timing models without fitting. We find that for even 
the best-timed pulsars to date, if standard model fitting 
is employed on data sets not significantly longer than 
1 year, the noise-induced errors in the astrometric param- 
eters can exceed the precisions achievable with VLBI in 
similar time spans. Furthermore, by incorporating VLBI 
astrometry into timing models for newly discovered pul- 
sars, after collecting data for just two years, the power 
present in post-fit residuals from red processes having 
spectral indices characteristic of a SB of GWs can be 
enhanced by a factor of 20. VLBI can additionally be 
used to aid multi- wavelength studies of the rapidly grow- 
ing population of MSPs; characterizing and monitoring 
these new MSPs will help to grow PTAs and thus aid in 
GW detection. 

For known pulsars with decades of timing data, VLBI 
astrometry has an advantageous but less significant influ- 
ence on the timing. It dramatically enhances the chances 
of detecting coherent GW signals with periods of six 
months or a year. Aside from this small class of GW 
sources, using VLBI astrometry will only boost the sig- 
nal of a 7 = 13/3 red process by about 2% in ten years of 
timing data. The reduction in the number of parameters 
being fitted in a timing model is a permanent improve- 
ment for tests for the presence of signal in the data, 
but this becomes less significant as more TOAs are col- 
lected and the number of degrees of freedom in the data 
grows. For some of the most well studied pulsars, those 
in relativistic binaries, even those having decades of tim- 
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ing data, VLBI astronictry will be of higher precision 
than timing astrometry and will allow correction for the 
Shklovskii effect to improve measurements of the PK pa- 
rameters of these systems. The decrease in the number of 
timing model parameters being fit for and the ability to 
augment the power of a SB of GWs in the timing residu- 
als of a pulsar in certain cases by a factor as substantial 
as 20 may play a key role in the initial direct detection 
by PTAs of GWs. 
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